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Abstract 

The NASA/GSFC Scanning Raman Lidar (SRL) was 
stationed on Andros Island in the Bahamas during Au- 
gust - September, 1 998 as a part of the third Convection 
and Moisture Experiment (CAMEX-3) which focussed 
on hurricane development and tracking. During the pe- 
riod August 21 - 24, hurricane Bonnie passed near An- 
dros Island and influenced the water vapor and cirrus 
cloud measurements acquired by the SRL. Two drying 
signatures related to the hurricane were recorded by the 
SRL and other sensors. Cirrus cloud optical depths (at 
351 nm) were also measured during this period. Op- 



deal depth values ranged from less than 0.01 to 1.5. 
The influence of multiple scattering on these optical 
depth measurements was studied. A correction tech- 
nique is presented which minimizes the influences of 
multiple scattering and derives information about cir- 
rus cloud optical and physical properties. The UV/IR 
cirrus cloud optical depth ratio was estimated based 
on a comparison of lidar and GOES measurements. 
Simple radiative transfer model calculations compared 
with GOES satellite brightness temperatures indicate 
that satellite radiances are significantly affected by the 
presence of cirrus clouds if IR optical depths are ap- 
proximately 0.005 or greater. Using the ISCCP de- 
tection threshold for cirrus clouds on the GOES data 
presented here, a high bias of up to 40% in the GOES 
precipitable water retrieval was found. 
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1 Introduction 


Raman lidar has long been regarded as one of the leading techniques for remotely quanti- 
fying numerous atmospheric parameters including water vapor, aerosols, temperature and 
clouds. Due to this broad measurement capability, the NASA/Goddard Space Flight Center 
(GSFC) Scanning Raman Lidar (SRL) was selected to participate in the NASA sponsored 
CAMEX-3 (third Convection and Moisture Experiment) hurricane study program which 
occurred during the months of August and September, 1998. The SRL was stationed on 
Andros Island, Bahamas during the experiment and acquired nearly daily measurements 
of water vapor, aerosols and clouds. SRL measurements of the variation of water vapor 
and cirrus clouds during the nearby passage of hurricane Bonnie from August 21 - 24 are 
presented here and constitute the first ground-based lidar water vapor and cirrus cloud mea- 
surements acquired in a hurricane environment. 

Significant drying episodes during the passage of hurricane Bonnie were observed and 
are likely due to mid-troposphere subsidence. The influence of multiple scattering on 
hurricane-induced cirrus cloud optical depth measurements was studied. A new cirrus 
cloud analysis technique will be presented which corrects for the influence of multiple scat- 
tering and also determines important optical and physical properties of the cirrus clouds. 
Cirrus cloud optical depths measured in the ultraviolet region of the spectrum are then 
translated to optical depths at the 1 1- and 12- micron channel location of the GOES satel- 
lite. Using these IR optical depth values, the influence of hurricane-induced cirrus clouds 

on GOES retrievals will be studied by comparing radiative transfer model simulations to re- 
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trieved surface temperatures and precipitable water. Using the International Satellite Cloud 
Climatology Project (ISCCP) [30] cirrus detection threshold on GOES data, the influence 
of undetected cirrus on GOES measurements will be studied. 

2 CAMEX-3 

Errors in prediction of hurricane track and thus landfall location are both dangerous for 
inhabited areas and can lead to unnecessary evacuation expense. And yet, small changes 
in initial atmospheric conditions can lead to large differences in the forecast of hurricane 
track and intensification [28], Because of this, CAMEX-3 was sponsored by NASA’s At- 
mospheric Dynamics and Remote Sensing Program with the goal of acquiring detailed 
measurements of water vapor, temperature and winds which can be used to help improve 
hurricane model initialization and forecasting. Several instrumented aircraft were sited 
at Patrick Air Force Base in Florida and made numerous flights in and near hurricanes 
Bonnie, Danielle, Earl and Georges as a part of this effort. Active and passive remote 
sensing instruments on board these aircraft were used to measure numerous atmospheric 
parameters including water vapor, winds, temperature, rainfall velocities and lightning 
(http://ghrc.msfc.nasa.gov/camex3/). 

A highly instrumented ground station was established on Andros Island in the Bahamas 

at the U. S. Navy’s Atlantic Undersea Test and Evaluation Center (AUTEC) as a part of 

CAMEX-3. Analysis of historical data indicated that the prevailing winds at AUTEC are 

out of the southeast during hurricane season. This indicated that this land-based loca- 
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tion on the windward coast of the island should give a good representation of the water 
vapor conditions over the open ocean. In addition to the SRL, this site included a Univer- 
sity of Wisconsin Atmospheric Emitted Radiance Interferometer (AERI) [14], radiosonde 
launches provided by both NASA/GSFC Wallops Flight Facility and University of Wiscon- 
sin, Global Positioning System (GPS) measurements of total precipitable water, Cimel sun 
photometer measurements of total precipitable water and aerosol optical depth at several 
wavelengths [29], chilled mirror hygrometer (http://www.humid.com/geiindex.html) mea- 
surements of relative humidity as well as standard ground measurements of temperature, 
pressure and relative humidity. 

The ground station served two main functions during CAMEX-3: 1) as a calibration 
and validation facility for CAMEX-3 and 2) as a source of highly detailed, long-term mea- 
surements of water vapor, aerosols, temperature and other parameters in the sub-tropics 
during hurricane season. Throughout the experiment, the research aircraft made numer- 
ous calibration/validation overflights of Andros Island allowing ground-based and airborne 
measurements of water vapor, temperature and winds to be compared. In this paper, we will 
describe a short segment of the nearly 2 months of measurements acquired at the ground 
facility: a 4-day sequence of water vapor and cirrus cloud measurements taken between 
August 21 and 24, 1998 when hurricane Bonnie was in the vicinity of Andros Island. We 
believe this combined set of measurements to be the highest quality water vapor and cirrus 
cloud data ever acquired in a hurricane environment. The Scanning Raman Lidar and the 

other water vapor measuring instruments used in this study will next be briefly described. 

5 



3 The Scanning Raman Lidar 


The Scanning Raman Lidar is a mobile lidar system designed to measure water vapor [23] 

[45] , aerosols [15] [17], cloud liquid water [24], cloud droplet radius and number density 

[46] , cloud base height [4] and upper tropospheric temperature [12], The SRL detects light 
backscattered by molecules and aerosols at the laser wavelength as well as Raman backscat- 
tered light from water vapor (3657 cm -1 ), nitrogen (2329 cm -1 ), and oxygen (1555 cm -1 ) 
molecules. The SRL employs two different lasers for its measurements; a XeF excimer 
laser (351 nm output) for optimized nighttime measurements and a tripled Nd:YAG laser 
(354.7 nm) for daytime measurements. The receiving telescope is a 0.76 m, F/5.2, vari- 
able field-of-view (0.25 - 2.5 milliradians) Dall-Kirkham system mounted horizontally on 
a 3.7m optical table. The telescope field-of-view is steered with a large (1.2m x 0.8m), 
motorized flat mirror which rotates on a horizontal axis and is also mounted on the optical 
table. The optical table can be slid out the back of the trailer to allow atmospheric profiles 
to be acquired at any angle in the plane perpendicular to the trailer or continuously scanned 
from horizon to horizon. Alternatively, the lidar system may be operated completely inside 
the trailer by directing the output laser beam through one of three windowed openings in 
the trailer. Use of these windows allows vertical measurements and measurements at 5-10 
degrees above the horizon in either direction to be acquired. It also allows measurements 
to be made during rainfall. All of the SRL instrumentation, including lasers, large aperture 
telescope and data acquisition electronics, is housed within a single environmentally con- 
trolled mobile trailer which also has separate areas for new experiment development and 
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work space for several experimenters to perform data acquisition and data analysis. More 
information on the lidar instrument has been published recently [46] and is available at 
our website http://virl.gsfc.nasa.gov/srl/index.htm. The SRL measurements acquired dur- 
ing the passage of Bonnie were made during the nighttime to maximize the signal-to-noise 
ratio of the data. 

3.1 SRL water vapor mixing ratio calibration 

While it is possible to calibrate a Raman lidar absolutely [34], our past calibration efforts 

have demonstrated that a careful selection of radiosonde data [16] along with the use of 

a nitrogen filter calibration transfer technique [40] [45], yields a stable lidar calibration 

constant. For a period of approximately 7 years, from the first field deployment of the SRL 

for the Spectral Radiance Experiment in Coffeyville, Kansas in 1991 [7] until the CAMEX- 

3 deployment in 1998, the calibration constant of the SRL determined by comparison with 

a selection of Vaisala radiosonde data varied only ±3%. The calibration constant is the 

number by which the ratio of water vapor and nitrogen lidar signals must be multiplied to 

obtain water vapor mixing ratio. Optical modifications were made to the SRL prior to the 

CAMEX-3 deployment which have changed the calibration constant. This fact, coupled 

with concerns about the calibration of the Vaisala radiosondes launched during CAMEX- 

3 [25], have necessitated a more careful examination of the SRL water vapor calibration. 

For the CAMEX-3 field campaign, we have implemented a new calibration technique [13], 

which assumes that the atmosphere is saturated at the base of a cloud. 
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3.1.1 Cloud base calibration technique 


Very frequently during the CAMEX-3 field campaign, small cumulus clouds developed at 
the top of the marine boundary layer that was present at the SRL site. The atmosphere 
below these clouds was typically well mixed and therefore was characterized by an ap- 
proximately constant water vapor mixing ratio. These facts permitted the SRL water vapor 
calibration to be derived using the SRL measurements of water vapor acquired just below 
the cloud. The saturation mixing ratio was calculated at cloud base using temperature and 
pressure from a simultaneously launched radiosonde. This saturation mixing ratio was then 
used to derive the calibration constant to convert the ratio of lidar signals into water vapor 
mixing ratio [45], 

On 23 separate occasions during the CAMEX-3 experiment, the SRL water vapor mea- 
surements were calibrated in this manner. The mean calibration constant calculated from 
these comparisons is approximately 12% higher than the value that had been used for the 
previous 7 years of experimentation. This new SRL calibration constant has been used to 
analyze the data presented here. 

For comparison, the SRL calibration constant was also determined in the traditional 

fashion using the Vaisala radiosondes launched during the experiment. The radiosonde 

measurements were first re-scaled to compensate for errors due to package contamination 

[25] [22], The mean Vaisala-derived and cloud base derived calibration constants agreed 

to much better than 1%, The standard deviation of the Vaisala-derived calibration constant 

was 5% while the standard deviation of the cloud base derived calibration constant was 
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3%. Therefore, during the CAMEX-3 field campaign, the new cloud base calibration tech- 
nique agreed well with the traditional radiosonde calibration technique and showed more 
consistent results. 

4 CAMEX-3 ground site water vapor instrumentation 

Total precipitable water vapor (TPW) measurements from several different instruments 
have been analyzed as a part of this study. These instruments are the SRL, Trimble SSi 
GPS (U. Wise ), Vaisala RS-80 radiosonde (U. Wise.), VIZ hygristor radiosonde (WFF), 
Cimel sun photometer (NASA/GSFC), GOES satellite, and a combined technique that uses 
the AERI (U. Wise.) and GOES. All of the ground-based instruments except the Cimel sun 
photometer were situated within a 1 00m radius approximately 1 km from the east coast 
of Andros Island. The sun photometer was located approximately 1 km west of the other 
instruments. The SRL measurements were limited to the nighttime periods while the sun 
photometer data were limited to daytime. The instruments use different techniques to make 
their measurements of TPW which can influence the values derived. The instruments and 
those techniques will be briefly summarized here. 

4.1 Atmospheric Emitted Radiance Interferometer (AERI) 

The University of Wisconsin AERI instrument [14] measures infrared radiation between 
approximately 3 and 20 microns with less than 1 wavenumber (cm -1 ) resolution using a 


Fourier transform infrared spectrometer. Radiance spectra acquired every 10 minutes are 
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transformed into vertical temperature and water vapor profiles by inverting the radiative 
transfer equation [36], The first guess water vapor solution is a hybrid profile using a 
statistical ensemble of radiosonde measurements for the boundary layer and the National 
Centers for Environmental Prediction (NCEP) GOES satellite profile above the boundary 
layer [39], The AERI retrievals that result are limited to an altitude of approximately 3 
km. To calculate TPW, the GOES water vapor profile is used above the height of the AERI 
retrieval [39], If coincident GOES retrievals are not available, the closest available retrieval 
is used. The AERI instrument that was deployed to Andros Island is similar to automated 
ones that have been installed at the Southern Great Plains (SGP) Site of the Department of 
Energy’s Atmospheric Radiation Measurements (ARM) Program [39], Based on an exten- 
sive comparison of AERI+GOES and the Microwave Radiometer at the DOE SGP site, the 
RMS difference between the total precipitable water retrievals from the two instruments 
was approximately 0.8 mm [33], 

4.2 Radiosondes 

VIZ (manufactured by Sippican, Inc) and Vaisala radiosondes were launched from the An- 
dros ground-site during CAMEX-3 and acquired profiles of relative humidity, temperature, 
pressure and winds. The VIZ water vapor sensor is a carbon hygristor which uses changes 
in resistance to determine relative humidity. This is the radiosonde that was the standard 
for the U. S. during the latter half of the 20th century. Data processing errors in retrieving 
relative humidities from these radiosondes have been discussed [41] and new algorithms 

addressing these limitations implemented [8], The Vaisala RS-80 radiosonde uses a thin 
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polymer film whose dielectric properties change as a function of the amount of water vapor. 
The changes in capacitance created by these changes in dielectric constant are converted 
into relative humidity. Since 1998, it has been the preferred radiosonde used by the U. 
S. weather service. Relative humidity errors due to packaging of the radiosondes have 
been found [25] and algorithms implemented to correct for package contamination [22], 
Radiosonde measurements of total precipitable water are calculated from the profile of 
relative humidity and have been characterized as being accurate to the 1 mm level [47], 

4.3 Cimel sun photometer 

The Cimel sun photometer is a solar tracking instrument that monitors direct and dif- 
fuse solar radiation from which aerosol optical thickness, aerosol size distribution, aerosol 
phase function and precipitable water vapor are retrieved [29], The sun photometer was 
deployed to Andros Island as a part of NASA’s AERONET (AErosol RObotic NETwork) 
effort (http://aeronet.gsfc. nasa.gov:8080/). The goal of this program is primarily to measure 
aerosol properties. Due to this, the precipitable water vapor measurements are currently be- 
lieved to have an error of approximately ±10% (Brent Holben and Tom Eck, NASA/GSFC, 
private communication, January, 2000). The total precipitable water retrievals from Cimel 
presented here use the standard processing algorithm based on Lowtran line strengths. Re- 
cent results by Giver et al [18] indicating errors in the Hitran-96 line strength database 
would reduce the Cimel retrievals of total precipitable water by 13% [32], Errors in line 
strengths such as described by Giver et al are problems for all optical instruments retrieving 

water vapor and is therefore actively being studied [32], 

11 



4.4 Global positioning system (GPS) 


The measurement of total precipitable water vapor using ground-based GPS receivers is 
accomplished by estimating the excess zenith-scaled signal delay caused by the neutral at- 
mosphere [47] [6], The measurement uses observations from all GPS satellites in view at 
a fixed site, and requires improved GPS satellite orbits and Earth orientation parameters 
that are supplied by any one of the International GPS Service Orbit Centers [3], The sig- 
nal delays are caused by changes in atmospheric refractivity associated with temperature, 
pressure, and water vapor along the paths of the signals within a radius of about 1 1 km of 
a site in the mid-latitudes. 

During CAMEX-3, two different software packages and improved GPS satellite orbits 
were used to estimate the zenith tropospheric delays from the data acquired at Andros 
Island. One estimate was made by the NOAA Forecast Systems Laboratory (FSL) in Boul- 
der, Colorado using GAMIT software developed by the Massachusetts Institute of Tech- 
nology and improved orbits provided by the Scripps Orbit and Permanent Array Center 
(SOPAC) at the Scripps Institution of Oceanography. Another estimate was made by the 
GPS Science and Technology (GST) program within the University Consortium for Atmo- 
spheric Research (UCAR) in Boulder, using Bernese software developed by the University 
of Bern and CODE Astronomical Institute orbits from the University of Bern. Using either 
technique, precipitable water values are determined at approximately 30 minute intervals. 
Comparisons of GPS derived TPW versus radiosonde have indicated mean differences of 

less than 1 mm with an RMS difference also of less than 1 mm [47], 
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5 Comparison of TPW measurements 


The relative calibration of these instruments (or in the case of GPS, the data processing 
techniques) has been studied for the Bonnie passage period of August 21 - 24. The mean 
differences in TPW during this 4-day sequence of data are shown in figure 1 where the 
GAMIT-processed GPS precipitable water vapor measurements were chosen as a baseline 
since they fell roughly in the middle of the distribution. The mean differences are calcu- 
lated using the number of comparisons shown in the figure legend. The error bars indicate 
the standard deviation of the differences. For all sensors except the radiosonde and the 
Cimel, 30 minute average data sets were used. The radiosonde produces profiles which 
take approximately 1 hour to acquire. The Cimel data frequency varied from a few minutes 
to more than 30 minutes therefore a strict half hour average was not always possible for the 
Cimel. All of the data shown later in figure 3 were used to determine the statistics for this 
plot. 

The Cimel sun photometer results are the wettest of the group with a high bias of ap- 
proximately 9% with respect to the GAMIT GPS. It should be noted, though, that this com- 
parison is based just on two days of measurements since the instrument was dismounted 
during the day on August 22 as a part of hurricane preparations at AUTEC. As discussed 
earlier, if the retrievals had been performed with the new line strengths according to Giver 
et al [18], the total precipitable water results should be 13% lower. This would change the 
Cimel results from a wet bias of 9% to a dry bias of 4%. 

Discounting the Cimel, the water vapor instruments’ relative calibrations agree to within 
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Figure 1 : Precipitable water vapor differences among the various sensors stationed at An- 

dros Island, Bahamas during the period August 21-24 when Hurricane Bonnie passed 
nearby. The arbitrarily chosen baseline for the comparison is the GAMIT-processed GPS 
data. The error bars plotted show the standard deviation of the differences with respect to 
the GAMIT baseline. 
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approximately +/- 3-4% or 1.5-2 mm. The SRL is approximately 1% wetter than the base- 
line GAMIT GPS data while the Bernese-processed GPS data were approximately 3% wet- 
ter than the baseline. The relative difference between the two methods of GPS processing 
is consistent with other investigations [27], It is thought to come from slight differences 
in data processing strategies. The VIZ and Vaisala radiosondes were within 1-2% of the 
baseline with the VIZ showing wetter measurements than the Vaisala. The AERI+GOES 
retrievals were the driest of the group with average values approximately 4% lower than 
the baseline. These AERI+GOES results seem to be consistent with others which indicated 
a tendency toward a dry bias as the TPW increases [33], 

The spread in the relative calibration of the instruments shown in figure 1 exceeds the 
claimed accuracy of many of the instruments. The results shown in the figure are an in- 
dication of the challenges inherent in accurate measurement of atmospheric precipitable 
water vapor. Based on previous studies of such differences [27] and discounting the Cimel 
measurements, this level of agreement is actually quite good, however. We are aware of no 
other long-term water vapor measurements of this quality acquired in the sub-tropics dur- 
ing the passage of a hurricane. While the uncertainties in these water vapor measurements 
presented here can translate into significant errors in radiative transfer calculations [7], 
they nonetheless represent a significant improvement over radiosondes alone for studying 
hurricane evolution. 

6 Precipitable water vapor measurements during the passage 

15 




Figure 2: GOES water vapor imagery at 0615 UT on the night of August 22, 1998. Andros 

Island is indicated by the white box. Notice the dry subsidence region to the west of the 
hurricane. 

of hurricane Bonnie 

Bonnie became a hurricane on the evening of August 22, 1998 at a point eastward of the 
Bahamian islands. The GOES water vapor image of Bonnie at 061 5 UT is shown in figure 
2. Andros Island is indicated by the white box. 

Over the next 4 days Bonnie followed a generally northwest track striking the mid- 

Atlantic coast of the United States on the evening of August 26. The point of closest 

approach of the center of the hurricane to Andros Island (24.7 N, 77.8 W) was at a distance 

of approximately 500 kilometers to the east northeast of Andros on the evening of August 

24. Figure 3 shows TPW measurements made by the ground-site instruments during the 

passage of hurricane Bonnie. As Bonnie approached Andros over the period of August 21 

- 22, there was a distinct drying indicated by all instruments during this period. Values 

16 


of precipitable water vapor changed from approximately 60 mm on August 21 to approx- 
imately 40 mm on August 22. We believe this drop in TPW to be due to compensating 
subsidence in the mid-troposphere due to the hurricane. During August 22 and most of Au- 
gust 23 the TPW shows a gradual moistening to values of approximately 50 mm by midday 
on the 23rd as the subsidence region moved to the west of Andros. Later on August 23 and 
into August 24, a developing wave disturbance over the gulf of Mexico blocked and re- 
versed the westward movement of the dry region. This resulted in a secondary dry feature 
in the Andros TPW measurements by 0000 UT on the 24th. 

While the general agreement in the TPW measurements reported by the various sensors 
is quite good, there are interesting discrepancies to mention. The high bias in the Cimel 
measurements is evident. Differences among the other instruments can be explained at least 
in part by the fact that the instruments use different techniques to make their measurements 
of precipitable water vapor which can influence the values derived. For example, in general 
the SRL and GPS values compare reasonably well. However, the SRL data acquired on 
August 22 are in general lower than either of the GPS retrievals. This may be due to the 
volume averaging that occurs as a result of using 6-10 more-or-less randomly distributed 
GPS satellites to measure the zenith-scaled tropospheric signal delay. Satellite imagery 
tends to support this conclusion as well. The GOES-8 images such as the one, shown in 
figure 2 indicate that for much of the night of August 22 UT, the Andros ground site was 
at the edge of the dry region. The SRL measurement of TPW was made directly over the 


ground site while the GPS averaged over a region which included more moist air from 

17 




21 22 23 24 25 

Time (days of August 1 998) 


Figure 3: Evolution of precipitable water vapor as measured by Raman lidar, 

AERI+GOES, GPS (with two different processing algorithms), two types of radiosonde 
and sun photometer. Two drying periods associated with mid-tropospheric subsidence are 
evident: early on August 22 and August 24 (UT). The sensors agree in the general trends 
but there are specific differences that can be attributed to measurement techniques. See text 
for details. 
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surrounding regions. 


Another example of a discrepancy between the GPS and the SRL that can be explained 
by the averaging volume used is seen in the data of August 24. Here a difference of up to 8 
mm is seen between the SRL and GPS TPW values. This can be explained by the presence 
of localized clouds and showers over the Andros ground site during the night of August 
24 which greatly increased the TPW measured by the lidar but which did not significantly 
affect the GPS measurements. 

A final interesting point to mention concerning this figure is the different structure re- 
vealed by the two GPS retrievals which are reported at approximately the same 30 minute 
intervals. This is due to different constraints used in the retrievals that determine how much 
the water vapor content can change in a short period of time. The GAMIT-processed re- 
trievals are much less constrained than the Bernese results, which allows for more structure 
in the GAM IT retrievals. 

6. 1 Water vapor evolution as a function of height 

The profile measurements of water vapor made by the SRL can be used to study the height 
dependence of the changes in precipitable water seen in figure 3. The SRL water vapor 
measurements have been divided into layers and integrated to yield the precipitable water 
vapor by layer. These results are shown in figure 4. The layers used are 0-1 km, 1-2 km, 
2-3 km, 3-4 km, 4-5 km and 5-8 km. 

In general, figure 4 shows that the 0-1 km layer changes very little during the 4-day 
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Figure 4: The precipitable water vapor measured by the lidar has been separated into the 

contributions due to the layers 0-1 km, 1-2 km, 2-3 km, 3-4 km, 4-5 km and 5-8 km. 
This quantifies the change in column water vapor as a function of height. There was little 
change in the 0-1 km layer throughout this 4 day period indicative of a marine boundary 
layer. Significant changes are evident, however, in other layers. 
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sequence while precipitable water vapor (PW) contributions from the other layers vary 
appreciably. This indicates that most of the boundary layer moisture was under local control 
due to evaporation from the ocean, which is characteristic of a marine boundary layer. By 
contrast, middle and upper tropospheric moisture was greatly influenced by the subsidence 
associated with hurricane Bonnie. For example, between the nights of August 21 and 22, 
the largest differences in precipitable water vapor occurred in the 2-3 km layer with values 
changing from 8-10 mm to 4-6 mm. An interesting exception to the depletion of PW at 
higher altitudes is seen in the 2-3 km layer on the night of August 24 at approximately 0400 
UT when rain influenced the local water vapor environment increasing PW values from 
approximately 5 mm to approximately 10 mm. This indicates that a significant amount of 
rain likely evaporated before striking the ground. (There is actually a small enhancement 
to the lidar PW measurements during rainfall due to Raman scattering from liquid in the 
rain droplets [4], We estimate that this effect increased the precipitable water vapor values 
in the rainfall by approximately 1-2 mm.) 

7 Cirrus cloud optical depth measurements 

Accurate measurements of sea surface temperature and total precipitable water vapor are 
needed to improve hurricane track and intensification forecasting. Satellites offer the best 
chance of providing operational data as input to hurricane models. However, it is well 
known that the presence of cirrus clouds can pose problems for satellite retrievals. This 

is because thin cirrus clouds, while having small infrared emissivities, can be very cold. 
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Emission from these clouds can cause significant changes in satellite radiances compared to 
a cloud-free scene. A comparison of Raman lidar cirrus cloud optical depth measurements 
with retrievals of surface temperature and total precipitable water vapor from GOES-8 will 
now be performed to study the influence of thin cirrus on these satellite retrievals. The 
same 4-day period associated with the passage of hurricane Bonnie will be considered. 

The technique for calculating cirrus cloud optical depth using Raman lidar will first 
be briefly described. Then the magnitude of the influence of multiple scattering on these 
calculations will be quantified using a retrieval technique that determines both optical and 
physical parameters of the cirrus clouds. The optical depth values obtained in the UV 
will then be translated to the IR. A simple radiative transfer model will then be used to 
quantify the anticipated radiance seen by GOES satellite under varying cirrus conditions. 
Comparisons of the predictions of this model with values derived using the split window 
technique [37] will then be presented. Lidar measured TPW will also be compared with 
TPW retrieved from GOES data. The International Satellite Cloud Climatology Project 
(ISCCP) [30] cloud screening technique will be applied to these GOES data to study the 
influence of undetected cirrus on GOES TPW retrievals. 

7.1 Optical depth assuming single scattering 

The optical depth calculation from Raman lidar is based on the molecular nitrogen (or 

oxygen) signal which shows enhanced attenuation due to the presence of a cirrus cloud. 

The amount of this attenuation can be converted to optical depth once the atmospheric 
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density is known. The single scattering equation which yields optical depth is obtained by 
integrating the equation for aerosol extinction [1] and can be written as 



[a(A L ,r) + a(X N ,r)} dr 


In 


Nn (r 2 ) -P(Ajy,ri) 
rl N n (n) P(X N ,r 2 ) 


rn 

Jr 1 


[ojjTioi ( X[j , 7*) T ( A/v j ^*)j dn 


( 1 ) 


where rq is below the cloud, r 2 is above the cloud, A/, is the laser wavelength (351.1 
nm), A;v is the wavelength of the Raman nitrogen signal (382.4 nm), a(A x ,r) is the cloud 
extinction coefficient as a function of wavelength and range, N n (r) is the number density 
of atmospheric nitrogen (using the full atmospheric number density is equivalent) as a func- 
tion of range, P (A at, r) is the lidar Raman nitrogen signal and a mo i(X x ,r) is the extinction 
coefficient due to molecular scattering obtained from radiosonde data. Lidar measurements 
in cloud-clear regions indicated that aerosols did not contribute to the optical depths mea- 
sured at cirrus altitudes. Also, at these wavelengths, gaseous absorption is negligible and 
need not be included. 

Equation 1 yields the two-way optical depth which is the fundamental quantity measured 
by the Raman lidar. To convert this to a one-way optical depth, the wavelength scaling of 
cloud particle scattering must be considered. Assuming an Angstrom coefficient of k = 0 
in the following equation 


a(^L,r) _ 
a(X N ,r) 
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which is valid for cirrus particles that are typically very large with respect to the laser 
wavelength of 35 1 nm, the one-way optical depth at 35 1 nm is just one half of the two way 
optical depth shown in equation 1. 

However, equation 1 does not account for any multiple scattering that may occur in the 
cloud. The influence of multiple scattering is mainly due to one or more forward scatter- 
ing events accompanied by a single backscatter event [10], Multiple scattering is much 
more likely when large particles are encountered because of the intense forward scattering 
diffraction peak associated with these particles. This forward scattered component is added 
back into the beam and decreases the apparent attenuation of the beam. Thus, the influence 
of multiple scattering is to decrease the optical depth measured by lidar compared to the 
actual value. Lidar parameters such as the telescope field of view (2 milliradians for the 
SRL) and the laser divergence (1 milliradian) also influence the multiple scattering compo- 
nent of the signal. The influence of multiple scattering on the Raman lidar measurements 
of optical depth during the hurricane Bonnie passage period of August 21-24 will now be 
studied. 

7.2 Multiple scattering calculations 

As mentioned above, multiple scattering is much more likely when large particles are en- 
countered because of the intense forward scattering diffraction peak associated with these 
particles. As the particle size increases, forward-scattered light is confined to an increas- 
ingly narrow angular cone. This makes it more likely that a photon that is scattered forward 
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in a first scattering event will interact with another particle (the second scattering event) and 
be backscattered within the field of view of the lidar receiver. Equation 1 was formulated 
for single scattering only where the assumption is made that, if a scattering event occurs, 
the photon is lost from the forward-propagating laser beam. Therefore, in the case of large 
particles, which can scatter a large number of photons in the direction of the laser beam, 
the use of the single scattering equations can lead to errors in the calculated quantities. 

Most of the quantities derived from Raman lidar data are based on ratios of lidar signals. 
In case of ratio measurements, multiple scattering influences the numerator and denomi- 
nator nearly equally and thus tends to cancel [42], Examples of these quantities are the 
water vapor mixing ratio, liquid water mixing ratio, aerosol scattering ratio and the aerosol 
backscatter coefficient. However, optical depth is calculated using only a single lidar signal 
(e.g. Raman nitrogen) and, in the case of large particles, can be significantly influenced by 
multiple scattering. 

7.2.1 Multiple scattering equations 

The influence of multiple scattering on lidar signals is related to the optical depth of the 
scattering medium, the size of particles that are doing the scattering, the range to the scat- 
tering volume and specific parameters of the lidar system in use. This can be seen in the 
formulation of the multiple scattering equations developed by Eloranta using a Gaussian 
approximation for the forward scattered diffraction peak [10], The ratio of double and triple 

scattering to single scattering can be expressed as 
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is the optical depth. In these equations, P n is the signal intensity due to nth order scattering, 
V n (tt, R ) /V\ (tt } R) is the ratio of phase functions in the backscatter direction for an nth 
- order scattered photon and a singly scattered photon. For Raman backscatter, this ratio is 
equal to 1.0 due to the broad nature of the molecular phase function near the backscatter 
direction. The telescope half angle field of view is p t , p t is the laser half angle divergence, 
d is the depth of penetration into the cloud determined by the location of the backscattering 
event, is the extinction coefficient, and 0 S is the 1/e diffraction peak angular half-width. 
These equations have been reformulated from the published versions [10] in a manner that 
allows for more efficient numerical calculation. The diffraction peak angular width for 
spheres can be calculated from the form of the scattering amplitude for a spherical aperture 


calculated from diffraction theory given by [2] 
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Here 0 is the scattering angle measured with respect to the forward direction, x is the 
size parameter of the spherical particle defined as the circumference divided by the radius 
and Ji is the first Bessel function of the first kind. The intensity of scattering versus angle 
is therefore given by 


/(0) = |5(0)| 2 (7) 

7.2.2 Cirrus cloud multiple scattering corrections 

Cirrus clouds were above the SRL site on Andros Island for much of the night on August 

23, 1998. There were characterized by quite cold temperatures ranging from -45C to -75C 

based on radiosonde measurements. A thick portion of the cirrus cloud was used to study 

the influence of multiple scattering on Raman lidar measurements of cirrus optical depth. 

The results are shown in figure 5. In this portion of the cirrus cloud, the measured optical 

depth (before correction for multiple scattering) was approximately 0.6. Both second and 

third order multiple scattering were calculated assuming constant particle radii of 5 microns 

and 20 microns. These particle dimensions were chosen based on the retrieved particle sizes 

that will be presented later. For these calculations, the cirrus optical depth was obtained 

from the lidar-derived cirrus backscatter coefficient, which is essentially uninfluenced by 

multiple scattering [42], using the following equation. 
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In equation 8, S is the ’’bulk” extinction/backscatter ratio between rq and r 2 in units of sr 
and (3 ( x ) is the cloud backscatter coefficient ( km~ 1 sr~ 1 ). In this context, the term ’’bulk” 
is used to refer to a mean value through a cloud layer. The results shown in figure 5 use a 
value of 20 for S. 

Cirrus cloud optical depth calculated using equation 1 require two reference altitudes, 
r i and r 2 . The first altitude, rq , must be below the cloud and the second, r 2 , must be above 
it. One of the interesting points to make about figure 5 is that the influence of multiple 
scattering on measurements of cirrus optical depth becomes smaller as the upper reference 
altitude is increased. This effect will be used here to correct measurements of optical depth 
for the influence of multiple scattering and determine additional parameters of the cirrus 
cloud such as bulk extinction to backscatter ratio and bulk particle radius. 

The influence of changes in the upper reference altitude on calculations of cirrus optical 
depth is demonstrated in figure 6 using data acquired during the night of August 23, 1998 
at Andros Island, Bahamas. Upper reference altitudes of 17 and 20 km have been used. 
The optical depth using r 2 = 20 km clearly indicates higher values which is consistent 
with a multiple scattering influence. Also plotted is the optical depth error (x 10) for the 17 
km calculation. The error in the optical depth calculation helps to explain why the optical 
depth calculations at 20 km sometimes are less than those at 17 km. 

The difference in optical depth calculated at 20 km and 17 km is shown in figure 7. 
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Figure 5: Multiple scattering calculations for a cirrus cloud measured on the night of Au- 

gust 23, 1998. Cirrus particles of 5 microns and 20 microns were simulated. An extinction 
to backscatter ratio of 20 was used. The ratio of nth order scattering to first order scattering 
is plotted along with the cloud backscatter coefficient. The backscatter coefficient has been 
multiplied by 10 for easier viewing. 
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Figure 6: Optical depth calculations for a cirrus cloud measured on the night of August 23, 

1998. Upper reference altitude (r 2 ) of 17 and 20 km have been used to demonstrate the 
influence of multiple scattering. As r 2 is increased, there is a general trend toward higher 
optical depths as expected. The error in the optical depth retrieval (xlO) is also plotted. 
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Figure 7: Optical depth difference using a 20 km upper reference altitude versus 17 km. 

In general, the 20 km calculations yield a higher optical depth which is consistent with 
multiple scattering, however the random error in the data can, on occasion, cause the 17 
km value to exceed the 20 km vlaue. 

In general, the results using r 2 = 20 km are higher than for r 2 = 17 km with the dif- 
ferences as large as 0.09. However at times around 3.0, 5.0, 8.0 and 10.0 UT, the 17 km 
extraction produces lower results. Nonetheless, these measurements on average provide a 
quantification of the multiple scattering influence that can be used to determine other cloud 
parameters. 
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7.2.3 Cirrus cloud retrievals 


The optical depth that is required in the multiple scattering calculations using equations 
3 and 4 is determined through the use of equation 8. An iterative technique has been 
developed to determine the correct S that uses the optical depth difference shown in figure 
7, equations 3 and 4, and an initial value of S that is determined by using the optical depth 
calculated with an upper reference altitude of 20 km. Each iteration computes a new value 
of S by correcting for the multiple scattering computed using the previous value of S. The 
algorithm converges very quickly since the original optical depth measurements (using 
r 2 = 20 km) are in error by typically less than 5% for these cirrus cloud measurements. 
The second iteration of the algorithm produces less than a 1% change in the value of the 
optical depth of the cloud and thus retrievals presented here use 2 iterations. The width 
of the forward-scattered component of the multiply-scattered light is also solved for in 
this technique. This is given by the 1/e half-width of the forward diffraction peak, in 
equations 3 and 4. With this value known, the radius of the sphere which possesses the 
same diffraction properties can be determined using equation 6. The results of this retrieval 
technique are shown in figure 8. 

The original uncorrected optical depths measured at 17 km and 20 km are shown in solid 
and dashed lines, respectively. These are the same data as in figure 6. The corrected optical 
depth resulting from the iterative technique is also shown. All three of these have been 
multiplied by 100 for display purposes. The retrieved values for bulk extinction/backscatter 


ratio have a mean value in this cloud of approximately 20. The bulk radius of the particle 
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Figure 8: The results of an iterative technique are shown. This technique uses the optical 

depths calculated at different altitudes and the lidar-derived cloud backscatter coefficient 
to simultaneously correct for the influence of multiple scattering on cloud optical depth 
as well as to determine the bulk extinction/backscatter ratio and the radius of the sphere 
with the same bulk diffraction properties as measured in the cloud. The uncorrected optical 
depth calculations performed with = 17 km (solid line) and 20 km (dash) shown in 
figure 6 are repeated here. The retrieved values of corrected optical depth, bulk extinction 
to backscatter ratio, and particle radius are also plotted. 
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that has the same diffraction width as observed in the cloud is also plotted and has a mean 
value of approximately 10 microns. 

There are several important points to make using this figure. The first is that this tech- 
nique demonstrates that Raman lidar measurements of cirrus cloud optical depth can be 
corrected for the influence of multiple scattering and that, for these cirrus clouds, this cor- 
rection is small when compared to the optical depth calculated at 20 km. The average 
correction to the 20 km optical depth is less than 5%. A second point is that the bulk 
extinction^ackscatter ratio of the cirrus cloud can be determined. The average value of ap- 
proximately 20 for this cirrus cloud is consistent with other lidar measurements that have 
been made using a technique similar to this [11] [9], It is interesting to note that ray-tracing 
calculations based on actual in-situ cirrus crystal measurements indicate a much broader 
range of extinction/backscatter ratios than has been measured by lidar [52] [11], The rea- 
son for this is that the lidar is sensitive to the light scattering properties of the crystals 
which are related more to the particle area than to the particle long dimension. The optical 
theorem can be used to clarify this. 

The optical theorem can be formulated as [20] 

= Y Im ' f ( k = k °^ ( 9 ) 

where a t is the total cross section including scattering and absorption, €q is the polarization 

state of the incoming photon, k is the scattered wavevector, k 0 is the wavevector scattered 

in the forward direction and f is the normalized amplitude of the scattered electromagnetic 
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field. 


It is clear from equation 9 that the amplitude of the forward-scattered diffraction peak 
is directly related to the total extinction of the particle. In fact, Babinet’s principle [20] 
[2], which is appropriate for large, absorbing particles, implies that an equal amount of 
incident energy is diffracted by the particle as is absorbed by the particle. This means 
that the total extinction cross section of such a particle is twice the geometrical area of 
the particle. This result agrees with Mie theory in the large particle limit. Therefore, the 
forward-scattered diffraction peak of the multiply scattered radiation, which determines the 
multiple scattering influence on the data and forms the basis of this particle size retrieval 
technique, is directly related to the cross sectional area of the particle doing the scattering. 

It should be noted that the cross sectional area of a ’’large” ice crystal can actually be 
quite small because one of the crystal axes will typically be very small. In the case, then, 
of randomly oriented crystals, the projected area of the crystals can be representative of 
particles much smaller than the long crystal dimension alone would imply. Other recent 
investigations have indicated that small dimensions may offer good representations of cer- 
tain cirrus crystal properties. Grenfell and Warren [19] have shown that calculations of 
cirrus multiple scattering can be estimated accurately by using a collection of equal-radius 
spheres to represent each cirrus crystal. The radius is chosen such that the ratio of volume 
to area (V/A) of the collection of spheres equals that of the crystal. The radius of the equal 
V/A sphere in the case of hexagonal columns is approximately equal to the radius of the 

short axis of the crystal. Thus, it is sensible that typical retrieved radii using the multiple 
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scattering technique described here should roughly correspond to the dimensions of the 
short crystal axis. In light of these considerations and the fact that small crystals are asso- 
ciated with cold cirrus clouds [51] [35] such as these, the small radii retrieved here seem 
reasonable. 

7.3 Hurricane Bonnie cirrus clouds 

Raman lidar measurements of cirrus cloud optical depth (at 35 1 nm) acquired at the Andros 
Island ground site for the nights of August 2 1 -24 are presented in figure 9. For the purposes 
of this figure, the values have been determined using 20 km as the upper reference altitude. 
Based on the analysis presented above, this approach reduces the influence of multiple 
scattering to a few percent at most. If higher accuracy results were desired, the full multiple 
scattering correction technique could be applied. A 10-minute running average of lidar data 
has been used for these calculations. The error bars plotted indicate the uncertainty of the 
measurement according to Poisson statistics. 

The lidar data were also used to calculate cirrus altitude and geometrical thickness. 
These measurements indicate cirrus cloud base height ranged from a minimum of approx- 
imately 9 km at 0530 UT on August 21 when the cloud thickness was approximately 5 km 
(optical depth approximately 1.5) to approximately 16 km at 0900 UT on August 24 when 
the thickness was less than 1 km (optical depth less than approximately 0.01). There were 
times when lidar optical depth measurements were not possible. These were due to rain 

(after 0600 on August 21), system filter changes (0200-0245 on August 21 and 0200-0330 
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on August 22), low clouds (0500 - 0600 on August 24). The other period of no data occurs 
after 0600 UT on August 22 and indicates that no cirrus were detected by the lidar during 
this time. 

Because of the range of optical depths covered, the measurements of August 23 provided 
a convenient dataset to test the sensitivity of satellite retrievals to the presence of cirrus 
clouds. On this night the measured optical depth at 351 nm ranged from a minimum of less 
than 0.01 to a maximum of approximately 0.7. (It is interesting to note that the optical depth 
limit (at 694 nm) for visual detection of cirrus during the daytime has been determined 
to be ^ 0.03 [31]). The lidar cloud backscatter coefficient image is shown in figure 10. 
Here the backscatter coefficient is shown using a log scale with values ranging between 
approximately 3 x 10~ 4 and 3 x 10 -2 ( km~ 1 sr ~ 1 ). 

7.3.1 Discussion of the radiative impact of cirrus clouds 

In order to estimate the radiative effects of these cirrus clouds, a simple radiative transfer 
model which accounts for surface emissivity, surface temperature, cloud emissivity and 
cloud temperature was used. The model equation is 

Raat = (1 — £c) P(^sat-Ts) + £ c (3 (A 3at ,T c ) (10) 

where R sat is the predicted satellite radiance (Wm~ 2 sr~ l fxm~ 1 ), e c is the cirrus cloud 

emissivity calculated from e c = (1 — e ~ Tc ) where r c is the cirrus infrared optical depth, 

e s is the surface emissivity, P(X,T) is the Planck function, A sat is the wavelength of the 

satellite instrument channel, T s is the surface radiating temperature and T c is the mean 
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Figure 9: Four night sequence of cirrus optical depth as measured by the Scanning Raman 

Lidar. Values reported are for 351 nm and have been calculated using 20 km as the upper 
reference altitude to minimize the influence of multiple scattering. The error bars reported 
are those calculated from Poisson statistics based on the strength of the lidar signal. Note 
that the optical depth scale changes for each of the plots. 
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Figure 10: SRL measurements of cirrus cloud backscatter coefficient on the night of Au- 

gust 23, 1998 at Andros Island. This image uses vertical pointing data only. 
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cirrus cloud radiating temperature. The first term in equation 10 is the surface contribution 
to the satellite radiance and the second term is the contribution due to the cirrus cloud. The 
satellite effective brightness temperature T sat is then obtained numerically from the Planck 
function for the value of R sa t • Averaging over the GOES 11- and 12-micron channel filter 
widths is required since the index of refraction of ice varies significantly in this region of 
the spectrum [44], 

The purpose of this equation is not to yield highly accurate values of satellite radiance 
but rather to study the influence of varying cirrus optical depths on those radiances. Thus 
the radiative contribution due to atmospheric TPW, which was roughly constant during 
the measurement period, was not included. For the model calculations of radiance using 
equation 10, the values used were: e s = 0.95, T s = 302 K obtained from GOES during a 
cloud-clear period, T c = 214 K obtained from radiosonde measurement. 

In order to use the lidar measured optical depths for infrared radiative transfer calcu- 
lations, the optical depths must be translated to the IR. In previous studies, the ratio of 
visible (532 nm from a Nd:YAG laser) to infrared cirrus optical depth has been shown to 
vary between approximately 1.6 and 2.4 [49] [5], This is an important ratio to quantify 
since it translate approximately into the shortwave/longwave forcing due to a cirrus cloud. 
The ratio depends on particle size and, due to the changing values of the index of refraction 
of ice, the exact spectral locations that are being compared. These studies have indicated 
that the values for 1 1 microns can be larger than for 12 microns. 

To study the ratio of UV/IR cirrus optical depths, the same approach described in Wylie 
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et. al. [49] was used where the cirrus IR optical depth r c can be approximated using the 
following equations. 


t c = - In (1 - e c 


£c = 


7^4 7^4 

^ clr _ sat 

T&r ~ T' 


In these equations, T clr is the blackbody brightness temperature for a clear GOES pixel. 
Using the GOES brightness temperatures during the night of August 23 which are shown 
later in figure 12, the ratio of UV (at 351 nm from the SRL) and IR (at the GOES 1 1 and 
12 micron channel positions) cirrus cloud optical depth was evaluated using equations 1 1 
and is shown in figure 1 1 . 


The mean values for the ratio of optical depths shown in figure 11 are 1.6 ± 0.6 at 1 1 
microns and 1 .4 ± 0.5. at 12 microns. While there is significant uncertainty in these values 
due to the small sample size, these results point to the conclusions that 1) the 1 1 micron 
ratio is larger than the 12 micron ratio which is consistent with the VIS/IR ratio studies 


mentioned earlier and 2) these UV/IR values are approximately 20% lower than the VIS/IR 
ratios. Both of these results can be related to the size of the ice particles in the cloud. 


Depending on the size of the crystal, ice particle extinction efficiency can change quite 

significantly between 1 1 and 12 microns due to the large changes in index of refraction of 

ice in this spectral region [51] [35], When small ice crystals are involved (10-20 microns 

in radius), the extinction efficiency at 12 microns is significandy larger than at 1 1 microns. 

This effect has been observed in thin cirrus clouds using Nimbus-4 [26] and airborne mea- 
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surements [35], Therefore, the results of figure 1 1, which show an 1 1 micron UV/IR optical 
depth ratio that is larger than the 12 micron ratio, can be an indication of small particle sizes 
in these cold clouds (-45C to -75C). For very large diameter crystals such as would be ex- 
pected in cirrus uncinus, the extinction efficiencies at 11 and 12 microns should be quite 
similar so that one would expect the two curves in figure 1 1 to overlay each other [35], It 
should be mentioned, however, that since smaller crystals are expected at the tops of the 
clouds, the conclusion of small ice crystals based on the IR data alone could be influenced 
by a top-of-cloud bias in the IR radiances. Because of this possible bias, the location of the 
instrument making the IR measurement of optical depth, whether on the ground [5], from 
aircraft [35] or from satellite [49], must be considered in the analysis. With this in mind, the 
difference between the UV/IR optical depths ratios determined here and the VIS/IR optical 
depth ratios determined before [5] [49] are likely related to both the particle sizes in the 
cirrus clouds that were studied and to the techniques used to derive the IR optical depths. 
We will study this ratio with more GOES comparisons in the future, but, for this study, 
we have used the values of 1.6 and 1.4 as the scaling factors to adjust the lidar measured 
optical depths to those appropriate for the GOES 1 1 and 12 micron channel positions. 

Figure 12 shows the brightness temperatures calculated from equation 10 and the Planck 
function using the parameters described above for both the 1 1 and 12 micron GOES chan- 
nels (long dash and dot-dot-dash lines, respectively). Also plotted are the actual GOES 1 1 
and 12 micron channel brightness temperatures (closed boxes and triangles). The slight 

high bias of the model results with respect to the GOES data is consistent with the atmo- 
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spheric contribution that was excluded from the model. The retrieved skin surface temper- 
ature using the split-window physical retrieval technique [37] are shown using open dia- 
monds. No cloud screening was performed in these retrievals. Therefore, the GOES bright- 
ness temperatures and the subsequent retrievals have the effects of cloud-contamination 
implicitly in them. This was done for the GOES pixel that contains the lidar location (on 
the left in the figure) and for the adjacent pixel 5 km to the east of this location (on the right) 
and thus completely over the ocean. (The cirrus optical depth is plotted for comparison in 
figure 13). 

There are several points that can be made from this figure. Despite the sampling is- 
sues relating to the comparison of 10-minute averages of lidar data and ~5 km satellite 
pixels, these simple model calculations capture the main features observed in the satellite 
brightness temperatures. Also, the pixels over land (left) and over water (right) show good 
general agreement indicating that the constant surface temperature assumption in the model 
retrievals is reasonable. That being the case, the third point is that, for both of these pix- 
els, the changing cirrus cloud optical depth is the dominant factor causing fluctuations in 
the satellite brightness temperatures. This influence lasts until approximately 1000 UT as 
indicated by the general slope in the model predictions toward higher brightness tempera- 
tures. Taking 1000 UT as an estimate of the first time during the measurement period when 
the satellite brightness temperatures were uninfluenced by the presence of cirrus, the IR 
optical depth threshold above which the presence of cirrus significantly influences GOES 

satellite brightness temperatures is estimated to be approximately 0.005-0.01 based on the 
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Figure 12: Comparison of GOES 1 1 and 12 micron channel brightness temperatures and 

model calculations for the satellite pixel directly over the lidar site (left) and the pixel ~5 km 
to the east of the site (right). Also plotted are the retrieved skin surface temperatures using 
the split-window technique. The model assumes constant surface and cloud temperatures. 
The two pixels show generally good correlation of features. It is evident that even very thin 
cirrus clouds influence satellite radiances. For comparison, the cirrus optical depth values 
are in the next figure. 
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ER-scaled lidar measurements. This value is approximately an order of magnitude lower 
than the cirrus optical depth detection threshold goal established within the EOS science 
plan [21] in the discussion of required satellite measurements. Thus, it seems apparent that 
the next generation of earth sensing instruments may have biases in their retrievals due to 
undetected cirrus. It is again interesting to note the -0.03 optical depth limit (at 694 nm) 
for visual detection [31], Using an approximate scaling factor of 2.0 [50], the IR optical 
depth limit for daytime visual cirrus cloud detection becomes -0.015. There is a suggestion 
in these results that even a cirrus cloud that can not be seen by the naked eye might still 
have a noticeable radiative impact on satellite measurements. 

The corresponding precipitable water retrievals using the split-window technique are 
shown in figure 13. In this figure, the lidar-derived precipitable water and the cirrus optical 
depth measurements (adjusted to the IR) are also plotted. The lidar measurements indi- 
cate that the TPW changed relatively little during the measurement period. All significant 
variation in the retrieved TPW from GOES is attributed to the influence of cirrus. 

It is clear from figure 13 that the cirrus-induced errors in the retrieval of TPW are larger 

than and in opposite direction to those in skin temperature shown in figure 12. Increases 

in cirrus optical depth depress the retrieved surface temperature and elevate the retrieved 

TPW. A simple explanation for this effect can be obtained by considering the adjustments 

in the derived values of surface temperature or precipitable water required to account for the 

change in radiance due to the presence of cirrus. Due to the T 4 dependence of blackbody 

radiant energy, small reductions in retrieved surface temperature can explain the reduced 
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Figure 13: Retrieved precipitable water from GOES for the satellite pixel directly over the 

lidar site (left) and the pixel 5 km to the east of the site. Also plotted are the SRL measure- 
ments of precipitable water and cirrus optical depth. The correlation between retrieved PW 
and cirrus optical depth is clear. 
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brightness temperatures of a cirrus cloud contaminated scene. However, large increases in 
precipi table water are required to bring about comparable reductions in brightness tempera- 
tures since the radiance due to the precipitable water is roughly linear with amount of water 
vapor. Also, the precipitable water is concentrated near the surface and is characterized by 
a mean radiating temperature that contrasts little with the surface temperature. 

It is interesting to compare these results with other cloud studies. In their study of VAS 
data, Wylie and Menzel [48] concluded that 50% of the clouds with IR optical depths of 
0.1 or less went undetected. In a more recent study based on HERS data [50], these same 
authors comment that the CO 2 slicing technique used on HIRS data allows the detection of 
cirrus for IR optical depths above approximately 0.05. 

In the most recent ISCCP cloud data products [30], cirrus detection was performed using 
a simple threshold technique based on brightness temperature. A cirrus cloud is indicated 
over land if the 1 1 fim brightness temperature is 4 A less than what is determined to be the 
cloud-clear value while over water this threshold is 2 K (D. Wylie, private communication, 
2000 ). 

Using the ISCCP detection threshold over water of 2 K with the GOES brightness tem- 
perature data presented in figure 12 yields the following results. The 2 K threshold corre- 
sponds to a cirrus cloud ER optical depth of approximately 0.05 in good agreement with the 
threshold determined from HIRS data. Referring to figure 13, the GOES PW retrieval (us- 
ing the right-hand pixel over the ocean) is approximately 20% larger than the SRL value for 

optical depths of 0.05. Over land, where the ISCCP cirrus detection threshold is 4 K, the 
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cirrus optical depth threshold and the error in TWP become approximately 0. 1 and 40%, 
respectively. This high bias is consistent with other investigations [47], If this example is 
representative of the current state of cloud detection algorithms, it seems evident that the 
probability that undetected cirrus are introducing a high bias into the precipitable water 
database is high. This high bias needs to be considered in the context of hurricane model 
predictions since this study indicates that thin cirrus are commonly found in the hurricane 
environment. 

8 Summary and Conclusions 

The NASA/GSFC Scanning Raman Lidar (SRL) was stationed at Andros Island, Bahamas 
during August - September, 1998 as a part of the third Convection and Moisture Experi- 
ment (CAMEX-3) hurricane study program. Lidar measurements of water vapor and cirrus 
clouds have been compared with various other sensors during the four day period of August 
21-24 when hurricane Bonnie passed near the island. The relative total precipitable water 
calibration of the instruments was compared where the SRL measurements were calibrated 
using a new cloud base calibration technique. The cloud base calibration value agrees very 
well with the calibration value derived from radiosonde but shows less random variation. 
The Cimel sun photometer was found to exhibit a wet bias of approximately 9% com- 
pared with the baseline measurement of precipitable water from the GPS using GAMIT 
processing. The general agreement of the TPW measurements of the other instruments was 

±3 — 4%. Differences between the GAMIT and Bernese processing of approximately 3% 
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with the Bernese technique producing wetter results are attributed to minor differences in 
data processing strategies. 

The evolution of the precipitable water vapor during this period was studied using all 
the water vapor sensors at Andros Island. The measurements revealed two drying episodes 
related to hurricane-induced mid-tropospheric subsidence. We believe these to be the first 
extended ground-based measurements of water vapor made in the near vicinity of a hurri- 
cane and a significant improvement over radiosonde measurements alone. Using the SRL 
profiling capability, the evolution of precipitable water was studied by layers indicating that 
the predominant changes in column water occurred above 1 km. The layer between the 
surface and 1 km was representative of a well-mixed manne boundaiy layer with relatively 
constant mixing ratio throughout the layer. There also was evidence of mid-tropospheric 
humidification due to rainfall on August 24. 

The evolution of cirrus cloud geometiy and optical depth were studied as well. The influ- 
ence of multiple scattering on the lidar measurements was studied. An iterative technique 
was presented which corrects for the influence of multiple scattering and allows cirrus cloud 
bulk extinction to backscatter ratio and particle radius to be determined. After converting 
the UV optical depths to IR optical depths based on a comparison of SRL optical depth and 
GOES brightness temperatures, the predictions of satellite brightness temperatures from a 
simple radiative transfer model were compared with actual GOES brightness temperatures. 
These predictions indicated that satellite radiances are noticeably affected for cirrus optical 

depths above approximately 0.005. Larger errors were induced in the retrieved precipitable 
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water than in the retrieved skin temperatures. Undetected cirrus should present a consis- 
tent high bias in GOES satellite retrievals of TPW. Using the cirrus cloud detection criteria 
of the most recent ISCCP analysis indicates this bias is up to 20% over water and 40% 
over land. Errors such as these could influence hurricane model initialization since cirrus 
clouds are abundant in the vicinity of a hurricane. Furthermore, cloud climatology studies 
based on SAGE II observations [43] have indicated frequencies of sub-visual cirrus near 
the tropical tropopause of up to 70%. This implies that the influence of undetected cirrus 
on satellite retrievals could be quite significant in tropical regions. 


An important conclusion of this effort is that satellite retrieval algorithms need to be able 
to detect the presence of cirrus clouds with IR optical depths as small as 0.005 in order to 
avoid significant influences on satellite radiances and thus potential errors in retrievals. This 
is an order of magnitude lower than the cirrus optical depth detection goal established in the 
EOS science plan. Improved satellite measurement strategies such as the 1.375 /im cirrus 
channel of the MODIS instrument (http:// modarch.gsfc.nasa.gov /MODIS/) on the recently 
launched and upcoming Terra and Aqua satellites are needed to improve satellite sensitivity 
to cirrus. However, the 1.375 fim channel is only effective at detecting cirrus during the 
daytime. Therefore, studies similar to that performed here are needed to determine the 
effectiveness of cirrus detection from satellite during both the daytime and the nighttime 
to determine if there are diurnal biases in the satellite precipitable water record due to 


undetected cirrus. 
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